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Magnetic fields in the absence of spiral density waves - 

NGC 4414 



(N 
O 

o 

(N 

< 

00 



> 
o 

00 

oo 

o 

o 

S3 

Oh: 

6 

CO 

c3 



M. Soida 1 , R. Beck 3 , M. Urbanik 1 , and J. Braine 2 

1 Astronomical Observatory, Jagiellonian University, ul. Orla 171, PL-30-244 Krakow, Poland 

2 Observatoire de Bordeaux, UMR 5804 CNRS/INSU, B.P. 89, F-33270 Floirac, France 

3 Max-Planck-Institut fur Radioastronomie, Auf dem Hiigel 69, D-53121 Bonn, Germany 

Received November 1st, 2001; accepted July 27th, 2002 

Abstract. We present three-frequency VLA observations of the flocculent spiral galaxy NGC 4414 made in order 
to study the magnetic field structure in absence of strong density wave flows. NGC 4414 shows a regular spiral 
pattern of observed polarization B-vectors with a radial component comparable in strength to the azimuthal 
one. The average pitch angle of the magnetic field is about 20°, similar to galaxies with a well-defined spiral 
pattern. This provides support for field generation by a turbulent dynamo without significant "contamination" 
from streaming motions in spiral arms. While the stellar light is very axisymmetric, the magnetic field structure 
shows a clear asymmetry with a stronger regular field and a smaller magnetic pitch angle in the northern disk. 
Extremely strong Faraday rotation is measured in the southern part of the disk, becoming Faraday thick at 6 cm. 
The distribution of Faraday rotation suggests a mixture of axisymmetric and higher-mode magnetic fields. The 
strong Faraday effects in the southern region suggest a much thicker magnetoionic disk and a higher content of 
diffuse ionized gas than in the northern disk portion. An elongation of the 20 cm total power emission is also seen 
towards the South. Although NGC 4414 is currently an isolated spiral, the asymmetries in the polarized radio 
emission may be sensitive tracers of previous encounters, including weak interactions which would chiefly affect 
the diffuse gas component without generating obvious long-term perturbations in the optical structure. 
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1. Introduction 

The relative role of small-scale velocity perturbations and 
galaxy-scale gas flows in determining the evolution and 
structure of galactic magnetic fields is one of the most 
hotly debated questions in studies of the evolution and 



Zweibel & Heiles 1997, Beck et al 



structure of galactic magnetic fields (e.g. Zweibel 1996, 

1996p . Large-scale 



galactic magnetic fields are known to display a coher- 
ent spiral-like pattern of polarization B-vectors (see Beck 
et al. 1996 for a review), indicating a substantial radial 
component capable of resisting the shear due to differ- 
ential rotation. The spiral magnetic field pattern could 
be produ ced b y the dyn amo mechanism (Wielebinski 
& Krause 1993 , Beck et al. 1996 ) in which the small-scale 
motions are constantly feeding the large-scale poloidal (i.e. 
radial and vertical) field. The radial field could also be pro- 
duced by continuous field stretching by large-scale flows 
due to density waves or bars (e.g. Otmianowska-Mazur 
& Chiba 1995). In choosing a flocculent spiral, where floc- 
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culent means without large-scale spiral structure, density 
waves and bars should play no significant role in deter- 
mining the magnetic field orientation. 

In well-studied nearby spiral galaxies it is extremely 
difficult to discern the effects of the turbulent dynamo 
from those due to processes in spiral arms. These objects 
have well-developed grand-design patterns with strong 
density wave compression and/or a high concentration 
of star formation in spiral arms. The first process may 
strongly modify the magnetic field by effects of compres- 
sion and gas flows along the arms (Otmianowska-Mazur 
& Chiba 1995). In some galaxies like M83 spiral- like com- 



pression regions (as traced by aligned dust filaments) 
fill the whole interarm space. Strong star formation in 
spi ral ar ms acts destructively on regular fields (Beck et 
al. |1996D . Differences in turbulent activity between the 
arms and the interarm region may give rise to strong con- 
centration of regular magnetic fields between the stellar 
arms, thereby strongly influencing the global field struc- 
ture (Rohde & Elstner p.998[ ). 



Knapik et al. ( 2000 ) first reported that flocculent 
galaxies may possess regular fields with a strong ra- 
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dial component. Their study was made at low resolution 
and only at one frequency. In this paper we present a 
high-resolution multifrequency polarization study of the 
flocculent galaxy NGC 4414 (Thornley & Mundy |1997| ). 
The velocity field shows no evidence for non-circular 
gas flows (Braine et al. 1993| , Sakamoto 1996, Thornley 



NGC 4414 8.44 GHz TP/PI B-vecI 

~1 I 



& Mundy p.997| ), NGC 4414 has a high surface gas density 



and is forming stars fairly intensively, much like m our 
own galaxy, but is in no way a starburst. All this ensures 
good conditions for building up a global magnetic field by 
turbulent processes in a way unaffected by large-scale gas 
flows and compressions. 

Table 1. Basic properties of NGC 4414 



Other names 

R. A. igso 

Decl.i95o 
Inclination 
Position Angle 
Distance 
Morphol. Type Sc 



PGC 40692 
UGC 7539 
12 h 23 m 57;8 
31°29'58'.'0 
55° 
155° 

19.2 Mpc 



Reference 

de Vaucouleurs et al. ( 
de Vaucouleurs et al. ( 
LEDA 
LEDA 



1991 
1991 



Thim (2000) 



LEDA 



2. Observations and data reduction 

The maps of total power and linearly polarized radio 
emission of NGC 4414 at three frequencies were ob- 
tained using the Very Large Array (VLA) of the National 
Radio Astronomy Observatory (NRAO)[} At 8.44 GHz 
and 4.86 GHz the compact D-array configuration has 
been used. At 1.415 GHz and 4.86 GHz observations were 
made with C-array. The on-source time was 1.3 hours and 
1.7 hours in the C-array at 1.415 GHz and 4.86 GHz, re- 
spectively. In the D-array the observing time was 14 hours 
at 8.44 GHz and 13 hours at 4.86 GHz. The receiver band- 
width was 50 MHz. 

The intensity scale was determined by observing 3C286 
and calibrating with the flux densities taken from Baars 
et al. ( 1977 ). The position angle of the linearly polar- 
ized emission was calibrated using the same source with 
an assumed position angle of 33°. The phase calibrator 
1219+285 was also used to determine telescope phases and 
the instrumental polarization. 

The data reduction has been performed using the 
AIPS data reduction package. The edited and calibrated 
visibility data were Fourier transformed to obtain maps 
of the Stokes parameters I, Q and U at all three fre- 
quencies. At 4.86 GHz the visibilities from both array 
configurations were combined. The maps were weighted 
with ROBUST=0 for the best resolution and sensitivity 
compromise, yielding synthesized half-power beam widths 




1 NRAO is a facility of National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 
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Fig. 1. Total power map of NGC 4414 at 8.44 GHz with 
polarization B- vec tors superimposed onto the Ha image 
from Pogge (1989 and priv. comm.). Uniform weighting 
yielding a resolution of 7" and r.m.s. noise of 9 /xJy/b.a. 
has been applied. Contour levels are 30, 70, 110, 150, 250, 
400, 550,. . .^Jy/b.a. 

(HPBW) of 7" and 11" at 8.44 GHz and 4.86 GHz, re- 
spectively. At 1.415 GHz our C-array map also weighted 
with ROBUST=0, has an original HPBW of 14". The nat- 
urally weighted maps, more sensitive to extended struc- 
tures, have HPBWs of 11" and 16" at 8.44 GHz and 
4.86 GHz, respectively. The Q and U maps were combined 
to get maps of the linearly polarized emission (corrected 
for the positive zero level offset) and of the position angle 
of apparent magnetic vectors (B-vectors). 

NGC 4414 has been also observed using the Effelsberg 
100-m MPIfR telescope to determine the total power flux 
at 2.695 GHz. Several small maps of the galaxy scanned 
either in R.A. or in declination were made, then the maps 
were averaged and the flux density was determined by 
integrating the final map in concentric rings. 

3. Results 

3.1. Total power emission 

The contour map of the total power brightness of 
NGC 4414 at 8.44 GHz is presented in Fig. || The to- 
tal power emission shows a ring-like distribution with 
a central depression roughly at the position of the hole 
in the CO and Ha emission (Thornley & Mundy 1997, 
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Fig. 2. The contour map of total power emission from 
NGC 4414 at 1.415 GHz. The contour levels are 2, 3, 4, 5, 
7, 9, 13, 20, 40, 60, 100, 150, 200 x 170/xJy/b.a. The map 
has been smoothed to a resolution of 20". The map is over- 
lai d upo n the distribution of CO(1-0) emission (Braine et 
al. 1993 ) shown as white contours on a grayscale back- 
ground 



Sakamoto 1996 , Braine et al. 1993| ). Our total power map 
also shows two maxima. The southern maximum coincides 
roughly with an optically bright clump and with a group of 
three large Hn regions (Pogge 1989). The peanut-shaped 
northern peak extends eastwards from the major axis and 
forms a bridge between two groups of H II regions. These 
brightness maxima have almost the same amplitude rel- 
ative to the rest of the disk in our total power maps at 
1.415 GHz and at 8.44 GHz. They are thus certainly non- 
thermal. In the outermost disk the radio emission in our 
total power maps at 8.44 and 4.86 GHz is asymmetric 
with a steeper brightness gradient along the western disk 
boundary (Fig. [j]). A small peak north of the centre coin- 
cides with a bright star-forming region and is not visible in 
our 1.415 GHz map so we suspect that it has an increased 
fraction of thermal emission. A faint radio ridge extends 
eastwards from the southern end of the major axis. It co- 
incides with weak optical and Ha extensions. At low fre- 
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Fig. 3. The total power map of NGC 4414 at 4.86 GHz 
with a resolution of 4" with a greyplot of the same quan- 
tity The contour levels are 3, 5, 7, 9, 11, 13, 15, 17, 19, 
21, 23, 25, 27 x 7 AiJy/b.a. 



quencies, the total power map at 1.415 GHz also shows a 
clear tail extending to the south (Fig. ||). 

Our C-array map of NGC 4414 at 4.86 GHz shows 
the brightest total power peaks with a resolution of 4" 
(Fig. ||). The total power hole in the disk centre is very 
deep at this resolution with no trace of any nuclear source. 
In the southern disk there is a good agreement between 
the total power peaks and bright Hn regions. In the outer 
parts of the NE disk a large Hn region at R. A. 1950 = 
12 h 23 m 57!9 Deci 95 o = +31°30'50" corresponds to a local 
radio emission maximum as well. However, in the bright 
peanut-shaped feature the association is less obvious: its 
radio-brightest part lies between nearby Ha peaks. 

The spectral index distribution in the disk of 
NGC 4414 is shown in Fig. f|. The flattening of the 
spectrum in the NE disk corresponds to a bright com- 
plex of Hn regions at R.A.i 950 = 12 h 23 m 57!9 Deci 950 = 
+31°30'50". The mean spectral index in the inner disk 
(within the radius of 1') is about a « 0.82 (S v oc v~ a ) 
to the North a « 0.79 to the South. The difference is in- 
significant (r.m.s variations are about 0.02) and may be 
due to the flattening to a « 0.76 at the position of large 
Hn regions in the southern disk. Thus, there is no clear 
evidence for a flatter spectrum of total power emission in 
the southern disk region. 
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Fig. 4. The spectral index distribution over the disk of 
NGC 4414 fitted to all three frequencies. The contour lev- 
els start from 0.55 with a step of 0.05. All the maps are 
convolved to a beam of 14" 



3.2. High-frequency polarized intensity 

The high resolution map of the polarized intensity at 
8.44 GHz is presented in Fig. ||. The brightest peak of 
polarized emission is located in the northern half of the 
disk, east of the bright ensemble of H II regions. An elon- 
gated polarized feature is also located in the southern disk, 
just outside of the optically bright arm-like patch. The 
peaks of high polarized intensity do not clearly avoid nor 
correspond to bright star-forming regions. Diffuse, more 
smoothly distributed polarized emission is visible as well. 

Despite the lack of organized optical or Ha spiral struc- 
ture the polarization B-vectors form a clear spiral pattern. 
Our 8.44 GHz polarization map at high resolution (Fig. ||) 
shows that the magnetic field orientation does not fluctu- 
ate from place to place, being coherent even at the scale 
of our smallest beam (650 pc). The global magnetic spi- 
ral pattern becomes fully visible when natural weighting 
is applied, lowering the resolution to 11" but allowing 
smooth, low-brightness regions to be seen (Fig. ||). The 
B-vectors are nearly azimuthal (pitch angle near 0°) close 
to the northern semi-axis but show a strong radial com- 
ponent (pitch angle up to 45°) in the southern disk. In 
the very inner disk the polarization B-vectors cross the 
galaxy's centre along the minor axis, forming an S-shaped 
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Fig. 5. Contour map of the polarized intensity of 
NGC 4414 at 8.44 GHz with a resolution of 7" (uni- 
form weighting for highest resolution) and r.m.s. noise of 
8 /iJy/b.a. with B-vectors proportional to the polarization 
degree, overlaid upon the optical image from HST (per- 
mission from Dr W. Freedman from Carnegie Institution). 
Contours are drawn every 25 /itJy/b.a. 



pattern. On each side of the centre polarization minima 
are visible, presumably caused by beam depolarization of 
strongly twisted magnetic fields, seen perpendicular in sky 
projection. 

The highest degrees of polarization (up to some 30% 
at 8.44 GHz) were found along the western disk boundary 
and in the faint eastern tail. The inner disk shows a signif- 
icant north-south asymmetry. The mean polarized bright- 
ness at 8.44 GHz and the mean polarization fraction are 
about 40% higher in the North than the South, where the 
Ha emission is slightly stronger. In particular, the Ha map 
(Pogge 1989 and priv. comm.) convolved to the beam of 
23" (the resolution of the CO (1-0) map) shows the bright- 
est peak in the southern disk. The mean Ha brightness, 
within the area delineated by the level of 30% of the peak 
value, is some 10% higher in the southern disk than in the 
northern one. In this respect NGC 4414 resembles nearby 
galaxies in which the highly polarized emission tends to 
avoid regions of high star formation (Beck et al. 1996). 
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Fig. 6. Contour map of the polarized intensity of 
NGC 4414 at 8.44 GHz with B-vectors proportional to 
the polarization degree of NGC 4414. Natural weighting 
yielding an r.m.s. noise of 7 /LtJy/b.a., a beam of 11" and a 
higher sensitivity to extended structures has been applied. 
Contours are plotted every 25 /LtJy/b.a. The vector of 10" 
length corresponds to the polarization degree of 50% 



3.3. Faraday rotation and depolarization 

NGC 4414 lies very close to the North Galactic Pole where 
the foreground Faraday rotation measure (RM) is smaller 
than 30 rad/m 2 (Simard-Normandin & Kronberg 1980Q . 
Using the polarized background source found in our maps 
at about 4' NE from NGC 4414 we determined the fore- 
ground rotation to be —4 ± 8 rad/m 2 . For this reason no 
correction for the foreground rotation was applied to our 
data. 

The changes of RM with azimuthal angle in the disk 
(Fig. 0) do not show clearly any simple singly-periodic or 
doubly-periodic variations traditionally attributed to ax- 
isymmetric (ASS) or bisymmetric (BSS) global field sym- 
1996j ). Instead, the azimuthal profiles 



metry (Beck et al. 
show a sudden jump at an azimuthal angle of about 200° 
with RM values in the inner disk changing abruptly from 
+600 rad/m 2 to —600 rad/m 2 . The rapid sign change oc- 
curs almost throughout the southern half of the disk. Such 
a jump usually means that the Faraday rotation angle 
exceeds 90°, corresponding at 4.86 GHz to RM= 7r/2A 2 
~ 600 rad/m 2 . It does not imply any reversal of the mag- 
netic field in NGC 4414 but signifies a "Faraday-thick 
regime" when the rotation angle exceeds 90° (Sokoloff 
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Fig. 7. Azimuthal profiles of the Faraday rotation mea- 
sures of NGC 4414 between 4.86 GHz and 8.44 GHz inte- 
grated in sectors of azimuthal width of 5° along rings 15" 
wide with the same inclination and position angle as the 
galaxy. The azimuthal angle runs counterclockwise from 
the northern end of major axis. The naturally weighted 
maps were convolved to a common beam of 16". Angles 
have been corrected to face-on 



et al. |1998|) . Such strong Faraday rotation is exceptional 
at 4.86 GHz and has not been observed in moderately 
inclined nearby galaxies which become "Faraday-thick" 
much below this frequency (Beck at al. 1996). The jump 
is most conspicuous at small galactocentric radii, becom- 
ing smoother and of smaller amplitude in the outer disk. 
The azimuthal profiles shown in Fig. [?] also show that RM 
values in the northern half of the disk are of order of 50-70 
rad/m 2 , much smaller than in the southern half, especially 
at the galactocentric radii < 45". 

Details of the distribution of Faraday rotation between 
8.44 GHz and 4.86 GHz are shown in Fig. |[ The rotation 
measure (RM) forms large intermixed domains of positive 
and negative values which do not correspond to any of the 
known global field symmetries. Generally, positive values 
dominate in the SW and western disk part with a small 
spot at the northern tip of the major axis. In large parts of 
the southern and SW disk the RM exceeds +100 rad/m 2 . 
A narrow region of very strong Faraday rotation, reaching 
+600 rad/m 2 is located in the southern disk close the de- 
polarized region south of the centre, jumping on its other 
side down to —600 rad/m 2 . The negative rotation mea- 
sures occupy the SE and eastern disk regions, extending 
to northwest. In the northern disk the rotation measures 
are small (< 100 rad/m 2 ) with interspersed domains of 
positive and negative sign. This region is also substan- 
tially polarized at 1.415 GHz, yielding similar values of 
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Fig. 8. The distribution of rotation measures of 
NGC 4414 between 4.86 GHz and 8.44 GHz, overlaid on 



the greyscale Ha image (Pogge 198G and priv. comm.) 
showing the brightest H n regions. Circles denote positive 
and crosses negative RM. The symbol sizes depend on the 
absolute value of RM as shown in the legend. The contours 
show negative (dashed) and positive (solid) values of RM 
with a step of 50 rad/m 2 . The naturally weighted maps 
at both frequencies were convolved to a common beam of 
16" 



RM determined between this frequency and 8.44 GHz. 
The above picture only weakly depends on the assump- 
tion concerning small (< 30 rad/m 2 ) foreground Faraday 
rotation. 

The map of Faraday depolarization (DP) between 8.44 
and 4.86 GHz (defined as the ratio of polarization de- 
gree between the lower and higher frequencies) is shown 
in Fig. ^. In the southern disk a region strongly depolar- 
ized at 4.86 GHz is present, with a depolarization factor 
DP < 0.2. It is associated with a jump in RM, as expected 
for a depolarization due to strong Faraday rotation inside 



the emitting region (sec Burn 1966, Sokoloff ct al. 1998) 
and extends to the SE disk region. The southern disk is 
completely depolarized at 1.415 GHz which confirms the 
strong Faraday effects. Another moderately depolarized 
region was found NW of the centre, close to the major 
axis. The DP there is about 0.45, rising to 0.5-0.6 when a 
correction for depolarization due to RM gradients across 
the beam is applied. On average, the DP between 4.86 
and 8.44 GHz in northern disk is about 0.7. This region 
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Fig. 9. The distribution of Faraday depolarization (DP) 
in NGC 4414 between 4.86 GHz and 8.44 GHz determined 
from naturally weighted maps convolved to 16". The con- 
tour levels are 0.2, 0.4, 0.6, 0.8, 1.0 



was even detected in polarization at 1.415 GHz thus, it is 
definitely Faraday-thin at 4.86 GHz. 

The regions of strong Faraday rotation or depolariza- 
tion show some association with the ionized g traced 
by the Ha line. Both the Faraday rotation jump and the 
depolarized area lie close to the complex of three large 
Hn regions and the jump south of the centre coincides 
with the strongest peak of Ha emission. The area depo- 
larized by more than 50% close to the northern major 
semi-axis lies close to two large Hn regions. However, 
a similar complex of ionized gas NW of the disk centre 
(around R.A.i 950 = 12 h 23 m 56!3 Deci 950 = +31°30T3") 
shows a rather low RM and is not strongly depolarized at 
4.86 GHz. 



4. Discussion 

4.1. Total magnetic field strength 

The flux densities of NGC 4414, compiled from data avail- 
able in the literature are shown in Tab. 2. The value at 
2695 MHz has been obtained by integrating our Effelsberg 
map in circular rings (the galaxy was only barely resolved) 
out to a radius of 6'. The flux density at 4.86 GHz was ob- 
tained from the integration of our VLA naturally weighted 
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Table 2. The integrated radio spectrum of NGC 4414 



Frequ- Flux 
ency density 
[GHz] [mJy] 



error References 



[mJy] 



0.408 
1.490 
2.695 
4.800 
4.850 
4.860 
10.55 
10.70 



560 
231 
134 
76 
78 
83 
58 
43 



± 25 
± 2 
± 30 



± 
± 
± 
± 
± 



Gioia & Grego rini (198(1 ) 
Condon at al. ( |l99C ) 
This work. 
Gioia at al. (|l982|) 



Condon at al. (|199l|) 

This work. 

Niklas at al. (|l995|) 
Gioia at al. (|1982|) 



total power map in elliptical rings with the inclination and 
position angle from Tab. 1 out to a radius of 3'. 

The best power-law fit to this data yields a— 0.77. 
Using this value and assuming energy equipartition, a 
lower limit to the cosmic-ray spectrum to be 300 MeV (see 
Beck 1991), a ratio of proton-to-electron energy density of 
100 (Pacholczyk |1970 ), and a nonthermal disk thickness 
of 2 kpc (full thickness corrected to face-on) , we compute 
the mean magnetic field strength within the radius de- 
lineated by 5% of the maximum brightness at 4.86 GHz 
to be < B t >= 15 ± 4 fiG. The uncertainty takes into 
account factor two variations in the proton-to-electron ra- 
tio, the disk thickness and the lower energy cutoff as well 
as, a variation in the thermal fraction at 4.86 GHz be- 
tween and 10% (allowed by spectral fits to the data in 
Tab. 2). The total magnetic field in NGC 4414 is compara- 
ble to that of activel y star -forming (non-starburst) spiral 
galaxies (Beck et al. 1996 ). Assuming the regular field to 
be parallel to the disk plane we obtain its mean strength 
< B reg >= 4±lfiG. The average ratio < B reg > / < B t > 
is 0.25. 



4.2. Distribution of the nonthermal emission 

Significant correlations between the gas density, star for- 
mation rate and magnetic fields were established by 
Bcrkhuijsen et al. (1993). In grand-design objects these 
quantities are affected by processes in spiral arms, like 
the accumulation of molecular gas in dust lanes followed 
by massive star formation and the magnetic field enhance- 
ment in compression regions. 

Azimuthal distributions of the 4.86 GHz total power 
brightness of NGC 4414 (showing a smaller contamination 
by thermal emission than that at 8.44 GHz and less af- 
fected by the zero-spacing problem than the C-array map 
at 1.415 GHz) are compared to the Ha and CO(1-0) emis- 
sion in Fig. [h]. At all radii there is a reasonable correspon- 
dence between the maxima of emission at 4.86 GHz and 
those in the Ha and CO lines. In the innermost disk region 
the correlation is less obvious and the maxima shifted, pos- 
sibly because of severe beam-smearing effects. At larger 
radii the azimuthal variations of the total power bright- 
ness follow well those in CO and in Ha. The maximum of 
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Fig. 10. Azimuthal variations of the total power bright- 
ness of NGC 4414 at 4.86 GHz, compared to those of 



CO(1-0) (Braine et al. |1993j ) and Ha (Pogge [1989]) aver- 
aged in sectors with 5° azimuthal width along concentric 
rings 15" wide. All the data were convolved to a com- 
mon beam of 23". The profiles were normalized to the 
mean value in each ring. The disk was rectified to face-on. 
Horizontal bars denote the beam size 

Thermal emission azimuthal profiles 




180 

azimuthal angle [deg] 



3 hi) 



Fig. 11. Azimuthal variations of the thermal emission 
of NGC 4414 at 8.44 GHz, averaged in sectors with 5° 
azimuthal width along concentric rings 15" wide. All the 
data were convolved to a common beam of 16". A non- 
thermal spectral index of 1 has been assumed. The disk 
was rectified to face-on. Thick horizontal bars denote the 
beam size 

the total power brightness at the azimuthal angle of about 
45° is visible both in CO and Ha. However, instead of a 
second maximum visible in CO and Ha at the azimuth 
of ~ 120° the radio emission shows only an asymmetric 
extension. It is possible that the CO and Ha minimum at 
the azimuth of 90° has been smoothed out by cosmic-ray 
propagation effects. All profiles exhibit a similar rise for 
azimuthal angles > 180°. However, the CO peak at an az- 
imuthal angle of about 340° in the rings outside a radius 
of 30" does not correspond to a similar feature in the Ha 
line. It has no counterpart in the radio domain so the Ha 
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emission is not caused by strong absorption of the optical 
radiation from hidden star-forming regions. Thus, the in- 
terrelations between gas density, star formation rate and 
magnetic fields also exist in flocculent galaxies, in which 
the distributions of gas, young stars and magnetic field 
are not strongly modulated by density waves. They have 
a global character, while small-scale details may show sub- 
stantial place-to-place differences. 

The central depression in the radio emission, gas den- 
sity and star formation deserves attention. The gas in this 
region is depleted during the early strong star formation 
necessary for the high stellar densities found in galaxy 
nuclei, fn density wave galaxies the gas is replenished by 
radial flows in spiral arms, which feed star- forming pro- 
cesses in the centre and inner disk. Some flocculent galax- 
ies show an a ccum ulation of CO gas in the central region 
(Regan et al. ^00l[ ) but this is not the case for NGC 4414 
(Sakamoto 1996 ). The central minimum of all species in 
NGC 4414 is an argument for the absence of dynamically 
important radial gas flows. 

4.3. Distribution of the thermal emission 

The interpretation of our Faraday rotation data needs in- 
formation on the distribution of thermal electron density 
independent of the Ha data which suffer from extinction 
effects. The azimuthal distribution of thermal emission 
from NGC 4414 at 8.44 GHz (Fig. p|) was determined 
from our spectral index map (Fig. W) assuming the non- 
thermal spectral slope a„t constant over the whole disk 
and equal to 1.0 - the mean observed value of a in the 
disk outskirts. The detailed value of a n t has little mean- 
ing for our discussion as it can only shift vertically the 
curves shown in Fig. |ll| without changing the shape of 
azimuthal variations of the thermal flux. 

The distribution of the thermal emission shows two 
maxima of similar amplitude at both ends of the major 
axis. The northern maximum at an azimuthal angle of 
about 30° coincides well with a similar peak of total ra- 
dio power, CO and Ha. The southern maximum is only 
slightly stronger and roughly corresponds to a second peak 
of all these species at an azimuth ~ 220°-240°. We do not 
find any strong excess of thermal emission in the south- 
ern disk which would naturally explain why the Faraday 
effects are so much stronger than in the northern disk. 

The lack of much higher thermal emission in the south- 
ern disk also means that this region does not exhibit strong 
star formation hidden by absorption in the Ha line. Thus 
we cannot explain the lower degree of polarization at high- 
frequency of the southern region by stronger magnetic field 
tangling generated by star-formation. The depolarization 
by differential Faraday rotation at 8.44 GHz reaches 0.68 
at the RM jump. However, such strong depolarization oc- 
cupies a small fraction of the disk and cannot account for 
the systematically lower mean polarization degree of the 
southern disk. 



4.4. The magnetic field structure 



The CO velocity field observed by Sakamoto ( 1996 ) 
does not show any perturbations similar to density 
wave streamings (Visser 1980), which would cause ra- 
dial stretching of magnetic field lines. As also stated by 
Thornley & Mundy (1997) the optical filaments or arm- 
lets which they detected in the near infrared do not cause 
strong departures from otherwise smooth HI and CO ve- 
locity fields. Though the armlets are likely to be dynamical 
features, they do not constitute a well-organized grand- 
design spiral pattern (Thornley & Mundy 1997). Small- 



scale phenomena like chaotic, turbulent gas motions may 
play an important role in the dynamics and star-forming 
processes in NGC 4414. 

Despite the lack of organized spiral arms, the B-vectors 
(rectified to face-on) (Fig. |l^) form a very regular spi- 
ral pattern. The magnetic pitch angle ip observed in 
NGC 4414 is about zero in the northern disk, increas- 
ing to about 40°-45° over the rest of the disk. This is also 
visible in the azimuth-ln(r) frame (Fig. |l^) in which a log- 
arithmic spiral would form a straight line inclined by the 
pitch angle. The infrared armlets detected by Thornley & 
Mundy (1997) show only limited similarity to the mag- 
netic field structure. They generally have smaller pitch 
angles in the northern disk and larger ones in the remain- 
ing disk portions. However, the magnetic field is aligned 
only with some portions of armlets while the polarization 
B-vectors apparently run across other armlet segments in 
all parts of the disk (Fig. |l2|, Fig. |l3|). For the deprojection 
to face-on, we assumed that both the armlets and mag- 
netic vectors lie in the galaxy plane. As demonstrated by 
Knapik et al. ( |2000D the presence of a three-dimensional 
magnetic field geometry as generated by the dynamo does 
not strongly affect the orientation of deprojected B-vectors 
even at inclinations as high as 60°. Thus, while the arm- 
lets may still have some influence upon the magnetic field 
structure, local processes like field amplification by the 
turbulent dynamo process need to be considered. 

The magnetic vectors crossing the disk centre are very 
suggestive of an admixture of non-axisymmetric magnetic 
field components, especially as the field running through 
the centre along the minor axis cannot be explained by the 
limited angular resolution (see Knapik et al. 2000). The 
hypothesis of a mixture of modes is supported by an analy- 
sis which simultaneously uses position angles of B-vectors 
at 8.44 GHz and 4.86 GHz (thus involving Faraday rota- 
tion effects), made kindly for us by Dr Andrew Fletcher 
(priv. comm). The azimuthal variations of magnetic pitch 
angles are best reproduced by a mixture of axisymmetric 
(the azimuthal magnetic wavenumber m=0), bisymmetric 
(m=l) and m=2 fields (Fig. [l3] middle panel). Because 
the magnetic field is a vector (in contrast to e.g. density), 
for the BSS configuration m=l is required to get opposite 
signs of magnetic field components at azimuthal angles 
differing by 180°. The amplitudes of the higher modes rel- 
ative to the axisymmetric field are 0.6 (m=l) and about 
0.3-0.5 (m=2, varying with the radius). The axisymmet- 
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Fig. 12. B-vectors of polarized in- 
tensity of NGC 4414 overlaid upon 
the optical image from HST (per- 
mission from Dr W. Freedman from 
Carnegie Institution), both recti- 
fied to face-on view. The optical 
image was digitally filtered to en- 
hance the spiral pattern. The spiral 
armlets discu ssed in Thornley and 
Mundy (1997) are marked by sym- 
bols 




Fig. 13. Top: the distribution of 
absolute values of pitch angles as 
a function of azimuthal angle and 
ln(r), r being the galactocentric dis- 
tance in arcmin. The azimuthal an- 
gle runs counterclockwise from the 
northern major semi-axis. Circles 
denote negative and stars positive 
values. The dashed line marks the 
radius beyond which NGC 4414 ro- 
tates differentially (Sakamoto 1996). 
Middle: the orientations of magnetic 
B-vectors for a mixture of m=0, 
1 and m—2 modes as described 
in the text, compared to the ob- 
served orientations of magnetic B- 
vectors in NGC 4414 at 8.44 GHz 
(bottom panel) . Crosses on the right 
side of the graph show the beam size 
for two different distances from the 
galaxy's centre. Solid lines show the 
spiral armlets discussed by Thornley 
& Mundy ([L997D 



ric field used in this analysis forms a trailing spiral with 
a pitch angle ip of about 22°, the bisymmetric one has 
ip of only 2° while the m=2 component has a large pitch 
angle whose absolute value varies between 40° and 80°. 
This mixture of modes well explains the highest polariza- 
tion extending from the northern major semi-axis along 
the western disk edge (see Fig. ^|) and the weakest polar- 
ized intensity in the SE disk. It also yields two asymmetric 
regions of a constant sign of magnetic field along the line 
of sight (responsible for Faraday rotation) with the sign 
reversal at the azimuthal angles of about 180° and 270°. 



This bears some resemblance to the RM sign distribution 
in NGC 4414 at radii < V (Fig. 111). 



As suggested by Chiba (|1993|) and Moss ( |1996a[ ) the 
generation of bisymmetric magnetic fields may be boosted 
by the parametric resonance in strong density waves. 
In this work we show the evidence for higher modes of 
the field structure without strong effects of spiral arms. 
The apparent contribution of an m=2 mode may indi- 
cate a modulation of the lower modes e.g. due to the 
non-axisymmetric distribution of thermal gas (see Fig. [ll]) 
and Faraday effects. Furthermore, the models discussed 
assume only a two-dimensional magnetic field parallel 
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Fig. 14. Distribution of the sign of the line-of-sight mag- 
netic field (responsible for Faraday effects) in our model 
(a) assuming a mixture of m=0, m=l and m=2 modes 
compared to (b) the sign of Faraday rotation between 
8.44 GHz and 4.86 GHz for NGC 4414. The Q and U 
data were convolved to a beam of 21" and restricted to 
signal stronger than 2a r.m.s. 

to the disk. Their further development requires a three- 
dimensional divergence-free magnetic field, computed us- 
ing elaborate MHD models. 

Though NGC 4414 is currently an isolated spiral (e.g. 
Braine et al. 



1997 ), the observed magnetic field asymme- 



tries may also signify external interactions. Asymmetric 
total power gradients, smaller magnetic and optical pitch 
angles and a higher polarization degree on one side of the 
disk were observ ed in the wind-swept spirals NGC 2276 
(Hummel & Beck |1995| ) and NGC 4254 (Soida et al. ^996 , 
Chyzy et al. in prep.). NGC 2276 also has a radio tail 



noted by Condon (1983). Pitch angle asymmetries and an 
enhancement of the polarization on one side of the disk by 
field compression are found in the tidally interacting spi- 
ral NGC 3627 (Soida et al. [2001 and references therein). 
We note also that the large-scale distribution of Hi pre- 
sented by Thornley & Mundy (1997) shows at interme- 
diate radii a slightly larger gas extent and slower density 
decrease with radius along the southern semi-axis, where 
wc find larger magnetic pitch angles. A similar coincidence 
is found in interacting spirals (Chyzy et al. in prep.). 

Isolated objects often show signs of past weak interac- 
tions, like retrograde stellar orbits or H I warps (e.g. Jore 
et al. 1996). Multiple "minor merging" effects may oc- 



cur during the galaxy lifetime (Haynes et al. 2000), some 



times producing large-scale asymmetries in spiral struc- 
ture (Zaritsky & Rix |1997 ) or short low-brightness stellar 
or mol ecula r tails (Aalto et al. 2001 ). Recently Kornreich 
et al. (|200 1|) found that while the optical morphology re- 
turns quickly to a symmetric appearance, the gas dynam- 
ics constitutes a long-lasting memory of past, weak per- 
turbations. 

NGC 4414 may be an example in which a distant in- 
teraction is remembered by its magnetic field. The inner 
disk is relatively unperturbed (see Sect. 4.2) and no cen- 



tral starburst nor a nuclear source are present. NGC 4414 
could have accreted diffuse gas, perturbing only its outer 
disk. As demonstrated using numerical models of galax- 
ies interacting with ambient gas clouds (Otmianowska- 
Mazur & Vollmer 2002), asymmetries in the polarized in- 
tensity, magnetic pitch angles and Hi may last for more 
than 1 Gyr. The magnetic field of NGC 4414 may thus 
"remember" a perturbation which could have occurred 
1 Gyr ago. We note finally that the asymmetries found 
in the Ha, Hi distribution and kinematics are very weak 
and only visible in low-surface brightness features. They 
are more obvious in the magnetic field properties like the 
pitch angle variations by 45° and the dramatic asymme- 
try in Faraday effects between the northern and southern 
disk. Multifrequency polarization observations may thus 
serve as a complementary tracer of weak perturbations, 
poorly visible in other domains. 

External influences could in principle explain the 
present-time radial magnetic field component. However, 
a single event cannot explain the observed strength of the 
regular magnetic fields (see Sect. 4.1) which require stable, 
continuous reproduction of their radial component over 
long timescales (e.g. Brandenburg & Urpin 1998). Without 
constant replenishment, the magnetic fields would diffuse 
out of the galaxy in a timescale of 10 8 yrs (Parker 1979). 
Although one can still imagine a quasi-continuous set of 
merging events which permanently regenerate the radial 
magnetic field, the turbulent dynamo seems a more rea- 
sonable alternative. The observations presented here are 
evidence for dynamo action with much less contamina- 
tion from spiral arms than in grand-design spirals, using 
a considerably better resolution than in our earlier work 
(Knapik et 



2001 ) 



To reproduce the strength of the Faraday rotation is 
by far a more complex issue. In the northern disk we got 
reasonable agreement between our estimate of the regu- 
lar magnetic field and standard assumptions concerning 
the ionized gas properties. We assume that about 80% of 
the observed thermal flux originates in H n regions mak- 
ing a small contribution to the Faraday rotation because 
of their small volume filling factor. Taking 20% of the 
thermal flux to be due to a diffuse component (Walterbos 
& Braun 1994 ) with a face-on full thickness of 700 pc and 
a volume filling factor of 0.04 (Fletcher et al. in prep.) we 
obtain for the equipartition regular field (assumed parallel 
to the disk plane) RM values of ~ 100 rad/m 2 , in reason- 
able agreement with observations. In the southern disk 
the apparent RM jump signifies a Faraday-thick regime 
and thus an intrinsic RM in excess of 500 rad/m 2 . We can 
reproduce this taking the observed thermal flux and as- 
suming the same regular magnetic field but assuming that 
50% (a factor 2.5 higher) of the observed thermal emission 
comes from the diffuse gas. We also require a filling factor 
of the diffuse ionized gas 2.5 times higher, equal to 0.1. 
Additionally we need a much thicker magnetoionic disk, 
with a face-on thickness of 3 kpc. 

The observed magnetic field structure in NGC 4414 
may be the result of the following processes: 
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The turbulent dynamo generates the global magnetic 
field with a substantial radial component and a mix- 
ture of modes. 

Past external interactions could enhance the regular 
field in the northern half of the disk, decreasing the 
magnetic pitch angles in this region and increasing 
them in the South, as found in galaxies known to be 



affected by external influences (Hummel & Beck 1995, 
Soida et al. 1996[ Chyzy et al. in prep.). The non- 
axisy mmetric fields may be boosted by interactions 
(Moss 1996b|) and/or from ram pressure e ffects as mod- 
eled by Otmianowska-Mazur & Vollmer ( |2002|) . 
- External interactions could thicken the magnetoionic 
disk and make it contain more diffuse ionized gas in 
the South than in the North. This could then make the 
southern disk of NGC 4414 Faraday-thick at 4.86 GHz. 
Such a phenomenon is so far unique among spiral 
galaxies. 

We have found that the magnetic field may be a very 
sensitive indicator of the behaviour of the diffuse compo- 
nent of the interstellar gas. In particular, it may trace signs 
of interactions which would pass unnoticed when observ- 
ing the stellar or Ha emission. Studying their magnetic 
fields together with deep mapping of their environment 
in H I, X-rays and in low-excitation spectral lines like C n, 
N II, O I or O in provides an important clue to their physics 
and evolution. 



5. Summary and conclusions 

We have performed a three-frequency VLA study of the 
flocculent galaxy NGC 4414 known to have extremely 
weak traces of optical spiral structure and no evidence for 
non-axisymmetric gas flows. The data were analyzed to- 
gether with CO (1-0) and Ha maps, yielding the following 
results: 

- The galaxy shows a bright total power disk with a cen- 
tral depression corresponding to the minimum in Ha 
and CO brightness, which may serve as an additional 
argument for the absence of radial gas flows. 

- The correlation of the nonthermal brightness with the 
CO and Ha emission also holds in a galaxy without 
strong spiral arms and thus it is not due to an accumu- 
lation of cold gas and star formation in density-wave 
compression regions. 

- Despite the lack of spiral arms and of non-azimuthal 
gas flows the galaxy shows a clear magnetic spiral pat- 
tern with a significant radial component, resembling 
that in grand-design galaxies. Dynamo action is the 
most likely source of a significant radial magnetic field 
component. Admixtures of bisymmetric field and pos- 
sibly higher modes are likely to exist even in a floccu- 
lent galaxy. 

- There is some correspondence between t he opt ical arm- 
lets discussed by Thornley & Mundy ( 1997 ) and lo- 
cal magnetic field orientations. The optical structure 



seems to follow the general global asymmetry of mag- 
netic pitch angles, while the magnetic lines seem to be 
aligned with only some segments of individual armlets. 
- NGC 4414 shows a strongly asymmetric distribution of 
Faraday rotation. Our observations are suggestive of a 
much thicker magnetoionic disk and an increased rela- 
tive content of diffuse ionized gas (compared to classi- 
cal Hn regions) in the southern disk. These asymme- 
tries, together with those in the magnetic pitch angles 
and total power asymmetries suggest some past inter- 
action (e.g. merging with dwarf galaxies or more likely 
accreting gas clouds) in the history of NGC 4414. 

In this work we show that the regular spiral magnetic 
pattern observed in all nearby galaxies does not require 
strong density wave action. We propose that the magnetic 
pattern in NGC 4414 is due to the dynamo process free 
from contamination by flows occurring in grand-design 
spiral arms. On the other hand, the picture of the mag- 
netic field in NGC 4414 is still far from being clear and no 
good description of the three-dimensional magnetic field 
geometry is yet available, except for some rough and very 
qualitative guesses. We suggest that observing the galac- 
tic magnetic field may provide a clue to the properties of 
rarefied ionized gas, contributing greatly to the Faraday 
rotation but little to the Ha emission. This gas phase may 
be more sensitive to external interactions than cither the 
molecular gas or H II regions associated with recent star 
formation. Thus, we propose to use the magnetic fields as 
an indicator of past external interactions working even in 
cases when the optical information does not indicate sig- 
nificant perturbations. We note finally that a more quan- 
titative magnetic field description in NGC 4414 may come 
from extensive modeling using more sophisticated multi- 
component field topologies and realistic distributions of 
thermal gas and relativistic electrons. A deep study of its 
environment in H I line and X-ray is highly desirable. Such 
a study is planned for the near future. 
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